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Abstract 

The search for one-dimensional electron systems with a giant Rashba-type spin splitting is of importance 
for the application of spin transport. Here we report, based on a first-principles density-functional theory 
calculation, that Bi zigzag chains formed on a heterogeneous GaAs(llO) surface have a giant spin splitting 
of surface states. This giant spin splitting is revealed to originate from spin-orbit-coupling (SOC) and 
electric dipole interaction that are significantly enhanced by (i) the asymmetric surface charge distribution 
due to the strong SOC-induced hybridization of the Bi px, Py, and p^ orbitals and (ii) the large out-of-plane 
and in-plane potential gradients generated by two geometrically and electronically inequivalent Bi atoms 
bonding to Ga and As atoms. The results demonstrate an important implication of the in-plane and out- 
of-plane asymmetry of the Bi/GaAs(110) interface system in producing the giant spin splitting with the 
in-plane and out-of-plane spin components. 
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I. INTRODUCTION 


Spintronics is one of the most growing researeh fields in eondensed matter physies^. Reeently, 
the Rashba-type spin splitting^, where the spin-orbit eoupling (SOC) lifts the spin degeneracy 
in the inversion-symmetry broken environments such as surfaces of solids or interfaces of het¬ 
erostructures, has drawn much attention because of the possibility of coherent spin manipulation 
without external magnetic field. After the first observation of such spin-split bands at the Au(l 11) 
surface^, a number of experimental and theoretical studies have intensively examined the Rashba- 
type spin splitting at various heavy-metal surfaces^, heavy-metal overlayers on surfaces^-—, and 
semiconductor heterostructures^. Especially, the search for one-dimensional (ID) spin-polarized 
electronic systems is desirable for manipulating spin carrier&i'i^. In this regard, ID nanowires 
formed on Si surfaces such as Au chains on stepped Si surfaces^, Pt nanowires on Si(llO)^, and 
Pb atom wires on Si(557)^ have been investigated to show a relatively larger Rashba-type spin 
splitting compared to two-dimensional (2D) electron systems. 

Although the Rashba-Bychkov model^ has been remarkably successful in explaining the phe¬ 
nomena of the Rashba spin splitting in terms of the effective magnetic filed experienced by elec¬ 
trons moving in the potential gradient perpendicular to the surface, it significantly underestimates 
the spin splitting as an energy scale of ~10~^ eV which is about 10^ times smaller— than exper¬ 
imental measurements^. In order to estimate quantitatively the size of the Rashba-type spin split¬ 
ting, several recent studies have taken into account other factors missed in the Rashba-Bychkov 
model, such as the atomic SOC effect^, the in-plane anisotropy of surface potential^i^, the 
asymmetry of surface charge distribution^""—, and the presence of orbital angular momentum at 
surfaces^-—. Since such various contributions to the Rashba-type spin splitting can be enhanced 
in ID electron systems compared to 2D electron systems^*^, it is interesting and challenging to 
search for the ID electron system with a giant Rashba-type spin splitting. Here, by choosing a 
prototypical ID electron system formed on a heterogeneous semiconductor substrate, we propose 
an ideal platform for exploration of the giant spin splitting. 

In this paper, we present a density-functional theory (DFT) study of the giant spin splitting in 
self-assembled Bi zigzag chains on GaAs(llO). This zigzag chain is composed of two geomet¬ 
rically and electronically inequivalent Bi atoms bonding to Ga and As atoms (see Fig. 
forming an interface at which the broken inversion symmetry produces the large potential gradi- 
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FIG. 1: (Color online) Perspective and top views of the optimized structure of Bi zigzag chains on a hetero¬ 
geneous GaAs(llO) surface. The dotted line indicates the unit cell. The a: and y axes point along the [001] 
and [110] directions, respectively. The large, medium, and small circles represent Bi, Ga, and As atoms, 
respectively. Two Bi atoms bonding to Ga and As atoms are labeled as Bioa and BIas, respectively. 

ents perpendicular and parallel to the surface. We find that the SOC enhances the hybridization 
of the Bi px, Py, and p^ orbitals as well as Ga and As p orbitals, giving rise to the asymmetry 
of surface charge distribution around the Bi atoms and the Ga and As substrate atoms up to the 
third deeper substrate layer. This significant asymmetric surface charge distribution together with 
the out-of-plane and in-plane surface electric fields increases the magnitudes of SOC and electric 
dipole interaction, leading to a giant spin splitting with the in-plane and out-of-plane spin compo¬ 
nents. The present findings will provide important implications for understanding the underlying 
driving forces behind the possible giant spin splitting in ID nanowires formed on heterogeneous 
semiconductor substrates. 


II. CALCULATIONAL METHOD 


The present DFT calculations were performed using the Vienna ab initio simulation package 
with the projector-augmented wave method and a plane wave basis set^^^. For the treatment of 
exchange-correlation energy, we employed the generalized-gradient approximation functional of 
Perdew-Burke-Emzerhof (PBE)2Z. The Bi/GaAs(110) substrate was modeled by a periodic slab 
geometry consisting of the seven Ga and As atomic layers with ~28 A of vacuum in between the 
slabs. The bottom of the GaAs substrate was passivated by pseudohydrogen atoms^^ with 0.75 or 
1.25 e. We employed a dipole correction that cancels the artificial electric field across the slab^^. 
The k-space integration was done with the 24x 36 Monkhorst-Pack meshes in the surface Brillouin 
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zones (SBZ) of the 1x1 unit cell. All atoms except the bottom two substrate layers were allowed 
to relax along the calculated forces until all the residual force components were less than 0.02 
eV/A. 


III. RESULTS 


We begin to optimize the atomic structure of the Bi/GaAs(l 10) surface using DFT calculation 
within the generalized gradient approximation of Perdew-Burke-Emzerhof (PBE) in the absence 
of SOC. The optimized structure is displayed in Eig. 1. The calculated bond lengths at the inter¬ 
face are = 3.00 A, ^Bi-Oa = 2.77 A, and Jbi-As = 2.77 A, in good agreement with previous 
theoretical^ and experimental data:^ (see Table IS of the Supplemental Material^). It is notice¬ 
able that two Bi atoms bonding to Ga and As atoms (designated as Bioa and BiAs in Eig. 1) are 
geometrically and electronically inequivalent because the Bi—Ga and Bi—As bonds have some 
different ionic characters due to a larger electronegativity of As compared to Ga and Bi atoms. 
Here, Bioa positions 0.11 A-higher than BIas, and it has a 5d core-level shift of 0.31 eV to lower 
binding energy relative to that of BiAs, in accordance with the measured value of ~0.4 eV from 
x-ray photoemission spectroscop)*^^*^. Due to these asymmetric characters of the Bi chains, the 
Bi/GaAs(110) interface produces the out-of-plane and in-plane potential gradients, as discussed 
below. Eigure 2(a) shows the calculated band structure of Bi/GaAs(l 10) with the band projection 
onto the Px, Py, and p^ orbitals of Bi. It is seen that the band dispersion of the highest occupied 
surface state (hereafter, designated as SS) is nearly flat along the EX and YM lines while broad 
along the FT and XM lines, indicating a ID electronic state^>^ along the Bi chains. As shown in 
Eig. 2(a) and 2(b), the orbital character of the SS state exhibits a strong k-dependence within the 
surface Brillouin zone (SBZ): e.g., the px or p^ character along the FT line, px along the FX line, 
Px, Py, or p^ along the XM line, and p^ along the MY line. As a consequence of such complex 
k-dependent orbital characters of the SS state, SOC easily induces a hybridization of the Bi Px, Py, 
and Pz orbitals, as demonstrated below. 

Next, we examine the effect of SOC on the geometry and band structure of Bi/GaAs(l 10) us¬ 
ing the PBE-i-SOC calculation. It is found that the inclusion of SOC changes <iBi Bi, '^Bi-Ga, and 
<^Bi-As by less than 0.04 A (see Table IS of the Supplemental Material^). Figure 2(c) shows the 
PBE-I-SOC band structure of Bi/GaAs(l 10)^. The spin degeneracy of the SS state as well as other 
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FIG. 2: (Color online) (a) Band structure and (b) Bi px, Py, and p^ orbital characters of the SS state obtained 
using PBE. The corresponding results obtained using PBE+SOC are given in (c) and (d). In (a) and (c), the 
bands projected onto Bi px, Py, and p^ orbitals are displayed. Here, the radii of circle are proportional to the 
weights of corresponding orbitals. The energy zero represents the Eermi level Ep. The inset in (a) shows 
the SBZ of the unit cell, while that in (c) magnifies the energy dispersion of the SS state along the X —)-M 
and X^—M directions. In (b) and (d), the px, Py, and p^ orbital components are mapped by using red, 
green, and blue color channels with their brightness, respectively. In (e), constant-energy contours around 
the X point, taken at 0.11 eV below Ep, are plotted with spin texture. Here, the arrows with increasing their 
brightness represent the SAMs with the negative, zero, and positive out-of-plane components, respectively. 

states is found to be lifted over the SBZ exeept at the high-symmetry points (i.e., F, X, M, and Y 
points). Obviously, the PBE-i-SOC band projeetion [Fig. 2(e)] and orbital eharaeter [Fig. 2(d)] of 
the SS state indicate a strong hybridization between the Bi px, Py, and p^ orbitals, compared to the 
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PBE cases [see Fig. 2(a) and 2(b)], leading to a eleetrie dipole interaetion that will be diseussed in 
detail later. The inset of Fig. 2(o) shows a elose up of the spin splitting near the X point along the 
XM direetion, illustrating the eharaeteristie dispersion of a Rashba-type spin splitting. For eom- 
parison of the size of the spin splitting with those of previously reported Rashba system9^-->^>^2s^^ 

we fit the ^-dependent dispersion of the spin-split subbands of the SS state along the kx or ky di- 

2 

reetion by using the Rashba Hamiltonian H = -i- f(pxz) ■S, where m* is the eleetron effeetive 

mass, Ur the Rashba parameter, and p and 5 are a momentum operator and a spin angular momen¬ 
tum (SAM) operator, respeetively. Here, the spin-dependent eigenvalues beeome £± = j^±aRk 
with the different values of m* and aR along the the kx and ky direetions. From the ealeulated band 
strueture, the eharaeteristie parameters of spin-split subbands sueh as the momentum offset XkR 
and the Rashba energy Er [see the inset of Fig. 2(o)] are extraeted to estimate m* and aR using 
the Rashba Hamiltonian. The present values of XkR, Er, and aR at the high-symmetry points are 
listed in Table I, together with those of previous Rashba systems^"^>^^>^. We find that aR = 4.94 
(2.27) eV A along A—j-M (F—J-T) is larger than 1.86 (1.09) eV A along A—J-F (F—;-A), indieating 
an enhaneed spin splitting along the y direetion parallel to Bi ehains^. Remarkably, the present 
values of aR are among the largest ones so far reported for Rashba spin-split systems (see Table I). 
We note that the eonstant energy eontours around the A point, taken at an energy of 0.11 eV below 
Ef, exhibit a pronouneed anisotropy between the AF and XM direetions, yielding eoneentrie el¬ 
lipses [see Fig. 2(e)]. Moreover, the eonstant energy eontours around the F point at energies lower 
than ~0.2 eV below Ef are open in the FA direetion (see Fig. IS of the Supplemental Material^) 
because the SS state has a flat band eharaeter along the FA line [see Fig. 2(e)]. 

On the experimental side, angle-resolved photoemission speetroseop)^ observed a nearly flat 
surfaee state along the FA line and the other surfaee state between the midpoint of the FA line and 
the X point. The former surfaee state is ~0.8 eV above the latter one. These dispersion features 
of the two surfaee states are generally similar to our band strueture obtained using the PBF or 
PBF-i-SOC ealeulation. In addition, aeeording to the inverse photoeleetron speetroseopy (IPS) 
experiment of MeFean and Himpsel^, the Bi monolayer on GaAs(l 10) produeed two unoeeupied 
surfaee states around 0.9 and 1.9 eV above the bulk valenee-band maximum (VBM) at the F point. 
On the other hand, the IPS experiment of Hu et al.— observed a pronouneed Bi-derived surfaee 
resonanee state around 1.25 eV above Ef at the F point. Noting that the band gap of the GaAs 
bulk is about 0.9 eV, the 1.25 eV-peak in the latter IPS data^ may eorrespond to the seeond peak 
in the former IPS data.— From our PBE (PBE-i-SOC) band strueture, it is found that the positions 
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TABLE I: Calculated Rashba parameters characterizing the spin splitting of the SS state around high sym¬ 
metry points of the SBZ, together with those of previously reported Rashba systems. The corresponding 
values obtained using HSE-i-SOG^ are also given in parentheses. 



AkR (A-i) 

Er (meV) 

as (eV A) 

Bi/GaAs( 110)-This 




r^x 

0.11(0.13) 

62(53) 

1.09(0.82) 

r^Y 

0.05(0.04) 

51(39) 

2.27(1.91) 

Y^r 

0.13(0.14) 

71(85) 

1.09(1.20) 

Y^M 

0.06(0.06) 

10(11) 

0.33(0.37) 

M^Y 

0.05(0.06) 

24(23) 

0.87(0.82) 

M^X 

0.06(0.08) 

19(25) 

0.66(0.63) 

X^M 

0.03(0.03) 

80(75) 

4.94(5.00) 

x^r 

0.06(0.06) 

60(31) 

1.86(1.09) 

Au/Si(557) (ID)^ 

0.05 

- 

- 

Pt/Si(110) (ID)^ 

0.12 

81 

1.36 

Bi/Si(lll) (1D)22 

0.17 

68 

0.80 

Au(lll) (2D)^ 

0.012 

2.1 

0.33 

Bi(lll) (2D)4 

0.05 

14 

0.55 

Bi/Ag(lll) (2D)^ 

0.13 

200 

3.05 

BiTel (3D)^ 

0.052 

100 

3.8 


of two lower unoccupied states at the F point are 1.12 (1.12) and 1.56 (1.45) eV above the bulk 
VBM, respectively, in reasonable agreement with previous IPS experiment a l ' ^ . 

Recently, it has become known that the asymmetric features of the surface states at surface 
atoms are crucial to determine the size of a Rashba spin splitting through SOC and electric dipole 
interaction. The former Hamiltonian is given by Hsoc = ^(VVxp) -5, where c is the velocity 
of light and V is a crystal potential, whereas the latter one represents the electrostatic energy of 
electric dipole moment in the surface electric field. In Fig. 3(a), we plot the xy planar-averaged 
electron charge densities (Pa and pp) of the spin-split SS state at ko = 0.03 [see the inset 
of Fig. 2(c)] away from the X point along the XM line. It is seen that Pa and pp exhibit some 
delocalization up to the third deeper GaAs substrate layer. This delocalization feature strikingly 


7 









FIG. 3: (Color online) (a) Planar-averaged electron charge densities for the SS state at [see the inset of Fig. 
2(c)], obtained using PBE and PBE+SOC. In (b), the planar-averaged crystal potential and ^^(Pa + Pp) 
(in Ryd atomic units) at ko are displayed along the z direction. 

contrasts with the spin-unpolarized case of Pa=p [see Fig. 3(a)] obtained using the PBE calculation 
without SOC, which shows a highly localized charge character around Bi atoms. Also we note that 
the PBE-i-SOC result for Pa and pp shows the broad maximum consisting of two peaks, possibly 
due to Bioa and BIas atoms. Thus, we can say that the SOC induces a large asymmetry of surface 
charge distribution, which in turn contributes to the Rashba spin splitting by AEsoc = 
f dz^^(pa + Pp)^. Eigure 3(b) displays not only the the xy planar-averaged crystal potential 
y as a function of z but also the calculated integrand of AEsoc for the spin-split eigenstates at 
ko- It is found that the magnitude of AEsoc can be dominant at the positions near Bi layer and 
the first GaAs substrate layer. On the other hand, the contribution of electric dipole interaction 
to the Rashba spin splitting can be expressed as AEd = —Ad ■ Eg, where Eg denotes the surface 
electric field and Ad is the dipole moment difference obtained from Pa and pp. Eor pa and Pp at 
ko, we obtain their difference Apa_p [see Eig. 3(a)] and then calculate Ad = —0.19 eA along the 
z direction. Since V around the Bi layer represents a highly deep, asymmetric quantum well [see 
Eig. 3(b)], a sizable magnitude of AEd can be expected by a very large surface electric field. It is 
noted that the asymmetric features of the heterogeneous GaAs substrate as well as the Bi chains 
composed of two geometrically and electronically inequivalent Bioa and BIas atoms produce an 
asymmetric in-plane surface charge distribution as well as an in-plane potential gradient along the 
X direction (see Eig. 2S of the Supplemental Material^), thereby contributing to the spin splitting 
through AEsoc and AEd. Thus, we can say that the SOC-induced asymmetries of the out-of- 
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plane and in-plane surfaee charge distributions together with the out-of-plane and in-plane surface 
electric fields lead to the giant spin splitting in the Bi/GaAs(l 10) surface system. 

Figure 2(e) shows the helical spin texture with the in-pane and out-of-plane spin components 
along the constant-energy contours around the X point, where the SAM direction rotates anti¬ 
clockwise (clockwise) at the inner (outer) contour. Note that the spins along the outer contour, 
located off from the XF or XM line, have the radial in-plane component (directing perpendicular 
to constant energy contour line), indicating that the spin splittings involve some contribution of 
the Dresselhaus effect^>^. Such Dresselhaus spin-orbit splittings may arise from the asymmetry 
of electrostatic potential or charge density in the Bi layer (composed of Bioa and BIas) and the 
heterogeneous GaAs substrate. In addition, the Bi/GaAs(110) system has one mirror-plane Oy 
symmetry with the xz plane (see Fig. 1), which is a combination of the proper rotation of 180° 
(about the y axis) with the inversion. Therefore, when the spin vectors move along the constant- 
energy contours from the irreducible part [i.e., the darkened area in the inset of Fig. 2(a)] to the 
neighboring one through the mirror-plane Oy symmetry, the Sx and components change their 
sign but the Sy component remains unchanged. Consequently, the spin vectors at the points cross¬ 
ing the rx line are oriented perpendicular to the mirror plane, and the whole spin texture along 
the constant-energy contours also satisfies the time-reversal symmetry that reverses simultaneously 
the wavevector and spin [see Fig. 2(e)]. 

It is noteworthy that the SOC-induced asymmetric surface charge distributions of the SS state 
originate from the admixture of the Bi px, Py, and p^ orbitals [see Fig. 2(c) and 2(d)] as well as their 
hybridization with Ga and As p orbitals (see Fig. 3S of the Supplemental Material). These orbital 
mixings in the SS state result in the formation of orbital angular momentum (0AM). Recently, Park 
et al.— pointed out that the existence of 0AM on the surfaces of high atomic number materials 
produces the electric dipole moment that interacts with the surface electric field, giving rise to 
a Rashba-type spin splitting. Indeed, as shown in Fig. 4(a), the 0AM is present in the upper 
subband of the SS state and is oriented mostly parallel to the corresponding SAM^, manifesting 
the interplay between orbital ordering and SOC for the Rashba-type spin splitting. It is interesting 
to note that, around the F point and near the FT line, the planar component of spin vectors is 
mostly composed of Sx- This feature of the spin texture is likely to be associated with the band 
structure of the SS state showing ID character along the y direction: i.e., the direction of SAM (or 
0AM) is perpendicular to the electron group velocity v = ^Vii£'(k) which nearly points in the y 
direction (see Fig. 4S of the Supplemental Material). Such a locking of SAM (or 0AM) to the 
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FIG. 4: (Color online) (a) SAM and 0AM textures of the upper subband of the SS state and (b) the total 
spin polarization l^l = ^Sx^ + Sy^ + and the component along the FF and XM lines. In (a), SAM is 
drawn with a larger arrow head than 0AM. Here, the SAM vectors with the negative, zero, and positive out- 
of-plane components are represented with increasing brightness of the arrow head. The arrows for 0AM 
represent only the in-plane components. The values in (b) are given in units of /r/2. 

crystal momentum that maximizes the magnitude of ean be utilized for the non-vanishing 

spin transport along the Bi ehains. Figure 4(b) shows the values of the total spin polarization |5| 
and the eomponent along the FT and XM lines. We find that along the XM line the magnitude 
of is comparable to those of the parallel eomponents. This sizable (or L^) component reflects 
the presence of the in-plane dipole moment and the in-plane potential gradient (see Fig 2S of the 
Supplemental Material^) generated by two eleetronieally different Bi atoms bonding to Ga and 
As atoms. 


IV, SUMMARY 


We have predieted for the first time that Bi zigzag ehains self-assembled on a heterogeneous 
GaAs(l 10) surfaee have a giant spin splitting with the in-plane and out-of-plane spin components. 
By means of the DFT caleulations, we revealed that this giant spin splitting originates from SOC 
and electrie dipole interaetion whieh are signifieantly enhanced by the large asymmetrie surfaee 
eharge distribution and the large out-of-plane and in-plane potential gradients. It was thus demon¬ 
strated that the in-plane and out-of-plane asymmetry present in the Bi/GaAs(l 10) interfaee system 
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plays an important role in the giant spin splitting. Our findings are antieipated to stimulate eurrent 
experimental and theoretical studies for exploration of the giant spin splitting in other ID electron 
systems formed on heterogeneous III-V semiconductor surfaces. 
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